Thermally Induced Local Failures in Quasi-One-Dimensional Systems: Collapse in 
Carbon Nanotubes, Necking in Nanowires and Opening of Bubbles in DNA. 



Cristiano Nisoli 1 , Douglas Abraham 1,2 , Turab Lookman 1 and Avadh Saxena 1 
1 Theoretical Division and Center for Nonlinear Studies, Los Alamos National Lab, Los Alamos NM 87545 USA 
2 Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road Oxford, 0X1 3NP England 

(Dated: October 20, 2009) 

We present a general framework to explore thermally activated failures in quasi one dimensional 
systems. We apply it to the collapse of carbon nanotubes, the formation of bottlenecks in nanowires, 
both of which limit conductance, and the opening of local regions or "bubbles" of base pairs in 
strands of DNA that are relevant for transcription and danaturation. We predict an exponential 
behavior for the probability of the opening of bubbles in DNA, the average distance between flat- 
tened regions of a nanotube or necking in a nanowire as a monotonically decreasing function of 
temperature, and compute a temperature below which these events become extremely rare. These 
findings are difficult to obtain numerically, however, they could be accessible experimentally. 

PACS numbers: 61.46.-w, 65.80.-g, 68.35.Rh, 87.14.gk 



Nanowires PQ, carbon nanotubes [5] and DNA fila- 
ments [5] are examples of Quasi One Dimensional (Q1D) 
mesoscopic systems characterized by one very long di- 
mension, on the molecular length scale, compared to the 
others that are of the order of nanometers. Their theoret- 
ical and technological interest is very high and motivates 
an analysis of the mechanisms governing their failure, 
yet numerical computations become impractical in the 
dilute regime. It has been experimentally observed and 
theoretically predicted that nanotubes of radii above a 
critical value of about 3 nm collapse globally along their 
length with a pancake-dumbell cross-section |U [5J O [7] ■ 
However, for carbon nanotubes of radii below the critical 
value, and therefore globally stable, we expect that ther- 
mal fluctuations will also cause collapse in local regions 
to impact their structural and transport properties. The 
necking of nanowires under thermal fluctuations has sim- 
ilar consequences. DNA filaments at physiological tem- 
perature are known to locally open to form the so called 
DNA bubbles [H |9], relevant to transcription and den- 
tauration processes. The size of these bubbles has been 
investigated by computationally intensive studies based 
on the Peyrard-Bishop-Dauxois model (PBD) [THIHI] at 
non-zero temperature |12j . however, a theoretical under- 
standing of their behavior is still lacking. 

We present here a very general statistical mechan- 
ics framework to study such thermally activated events 
which is applicable to a variety of Q1D systems, and we 
show how to implement it practically to gain some insight 
on the phenomena described above. We predict an expo- 
nential behavior for bubbles in DNA (which has recently 
emerged from numerical model simulations but has so 
far not been accessible experimentally) and the depen- 
dence of the average distance between collapsed regions 
in a nanotube or necking in a nanowire, which effectively 
infinite below a critical temperature and monotonically 
decreases to a non-zero minimum at very high tempera- 
tures. 



Following previous work |13j . we consider a Q1D sys- 
tem like the one depicted in Fig. [T] and, for simplicity, 
restrict ourselves to one degree of freedom, x(l), where 
I is the coordinate along the length of the wire. The 
boundaries are allowed to oscillate under an elastic en- 
ergy |x' 2 (Z) and they interact via a local potential V(x). 
For instance, for a collapsing carbon nanotube, x(l) rep- 
resents the shorter axis of the collapse, as we described 
the collapsed posrtion as a "pancake" of heigh x. In 
models for solitons in DNA, e.g. the celebrated Peyrard- 
Bishop model [10] . x represents the distance between 
complementary bases, whereas V(x) is a Morse poten- 
tial, chosen because exactly solvable, although a simple 
square well would suffice |14j . In strained, 2-dimensional 
nano-islands, V(x) represents the elastic interaction be- 
tween the edges. The locality of the potential is often an 
approximation which is valid when the typical length oc- 
cupied by the collapse failure or instability is longer than 
the lateral size of the wire. The problems we describe 
thus posses a characteristic length, l c below which col- 
lapse cannot happen, simply because there is not enough 
available space. In a carbon nanotube l c is the length 
needed to go from a cylindrical to a collapsed configura- 
tion and back. Different is the case of DNA: we will be 
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FIG. 1: A schematic depiction of localized failure in a Q1D 
system, where x c defines the collapse width, d is the lateral 
size at equilibrium, and x(l) is the relevant degree of freedom 
for our problem, allowed to fluctuate in the coordinate I, < 
I < L. 
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interested in studying the length of bubbles, and thus we 
will consider "defect" the joined double filament configu- 
ration: because of that l c = 3.3 A, the distance between 
bases, but there are no locality issues at play, and the 
potential is genuinely local, the chemical bond between 
complementary bases. In nanowires, x reppresent their 
lateral size. 

This approach provides us with a natural discretiza- 
tion length for our problem: that, as we shall see later, 
considerably simplifies matters. For specificity we define 
"failure" as x dropping below some critical value x Cl a 
value that is suggested by the problem itself and might 
have nothing to do with the particular shape of V(x). For 
a nanowire it will be the threshold below which band gaps 
open and conductivity drops. For a carbon nanotube, it 
might be the distance corresponding to the metastable 
flattened state. For bubble in DNA it will be the dis- 
tance at which bases are linked by hydrogen bonds. 

With these assumptions, the partition function for a 
piece of the wire of lenght nl ci starting and ending in a 
failure is given by a path integral over the allowed fluc- 
tuations 



Z n = / Vx(l) 
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(1) 



with a;(0) = x(n l c ) < x c , and n integer. The normal- 
ization in T>x(t) is chosen so that standard transforma- 
tions |16l 117] allow us to write Z n as density matrix 



Zn x c 



of the one-particle, temperature dependent hamiltonian 
given by 
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(3) 



Now, the probability that a system of length nl c will not 
fail anywhere in the middle is given by 



Pn = 



(4) 



where £l n is the sum over fluctuations such that x(l) > 
x c : that is, the sum over the fluctuations that only lead 
to failure at the boundaries, or 

Q n = [ Vx(l) e -/5/o"Mi* ,3 «+n*«)] , (5) 

J x(l)>x a 

again with x(0) = x(nl c ) < x c . 

We divide all the fluctuations in the definition of Z n 
into n groups: those that never lead to failure, those that 
lead to failure after the first n — 1, n — 2, . . . , 2, 1 steps. 
In this way Q, n enters the convolution equation 
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reminiscent of the work of M. E. Fisher et al. [15], The 
vertex term v, a length, is introduced to account for the 
path integration between and x c around each ml c . For 
small x c one expects v = Xx c . 

We can solve ^ for fi n by deconvolving via the gener- 
ating functions Z(u) = Yl^Li z nU n , f2(u) = J2n°=i 
thus obtaining 



ft(u) 
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(7) 



from which J7„ can be derived by power expansion of 
Cl(u). For instance, by direct differentiation of ([7]), and 
using Z(0) = one immediately finds 



Pi = I- 



(8) 



As expected by definition, a piece of length l c cannot fail. 

For simplicity, let us assume (justifications follow later) 
that only the temperature dependent effective Hamilto- 
nian in equation ^ has a bound ground state ?po(x) of 
energy e D (functions of temperature) and that Z n reduces 
to a projector on that ground state, Then we can write 
Z n = (1 — p)"Uq"-!; -1 where 



1-p: 



vx n 



(9) 



and Mo = exp(eo). Summing simple series and following 
the procedures described above one finds f2 n = w _1 (l — 



dx / dy(x I e- nH e I y) (2) p)p r ' 



and immediately 

Pn = P" 



(10) 



where, correctly, P\ = \ and thus again, a piece of length 
l c does not fail. This exponential behavior of (10) has 



been found in numerical studies of bubbles in DNA, as 
explained later [12 . From the average n one finds the 
average distance / = l c h between failures for a long wire 



lr 
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(11) 



which shows that the lower is 1 — p (related to the proba- 
bility of penetration below x c of the eigenstate) the longer 
is the distance between failures. Also the probability P n , 
and hence I or n are independent of uq and thus of the 
actual eigenvalue of the state, as expected, since the prob- 
lem is invariant under uniform shift of the spectrum of 
(J3j): of course when more than one eigenvalue is relevant, 
then only the differences among them plays a role. 

When can we apply the outmost bound eigenvalue ap- 
proximation? We leave to the reader to prove that from 
([7]), when the entire spectrum of ^ is retained, f)(u) is 
a meromorphic function with real positive poles of or- 
der one. In the approximation of infrequent failures, or 
large bubbles in DNA (acceptable at physiological tem- 
perature) where the average distance n is large, only the 
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FIG. 2: Top: Plot of V(x) in Q. The system is in equilibrium 
between d and d— A, with a hard core repulsion at d and zero, 
and a square potential barrier at d — A of strength e . Dashed 
lines are eigenstate of eigenvalue larger and smaller than fie},- 
Bottom: the behavior of 1 — p(T) oc l^ 1 , for the potential 
above: notice saturation at high temperature, whereas defects 
are virtually inexistent at for T/T c < .8. 



smallest pole is relevant. Then, if the spectrum of ^ is 
such that e s — eo 3> 1, one only considers the bound 
groundstate, as the correction to p is proportional to 
exp(eo — €i). We will see that in the cases we study 
both approximations are fulfilled. 

For specificity we consider now the steplike potential 
depicted in Fig. [2j designed to capture in a simplified 
picture the essential features of carbon nanotube col- 
lapse and the rich phenomenology of nanowire necking 
(whereas a different approach will be described later for 
DNA bubbles). We assume that the system is in equi- 
librium between distances d,d — A, with A <C d as the 
inaccessible region is much wider than the size of the 
Q1D system. The infinite repulsion at d ensures stability 
of the wire toward large fluctuations, and, in the case of 
carbon nanotubes, it accounts for the very large energy 
needed to break sp 2 bonds, whereas the finite height ej, 
of the square repulsion allows for some penetration and 
thus a non zero probability of localized failure 1 — p in ^ . 
Hard core repulsion at zero prevents negative values of 
x, which would lead to an artificial additional extensive 
entropy. 

Now, a look at (|3| shows that for very small T, the 
potential barrier (3eb is very large in the Schrodinger op- 



erator, and thus the probability 1 — p from Q is very 
small, because the ground state provides little tunneling. 
As the temperature increases, the barrier decreases, the 
ground state has larger probability for penetration and 
thus 1 — p also increases. There is a definite tempera- 
ture T c such that when T > T c , there are no more bound 
states of eigenvalue less than (3tb and the behavior of the 
ground state in the forbidden zone becomes sinusoidal 
rather than exponential, and thus 1 — p becomes much 
larger. Yet 1 — p does not increase indefinitely with tem- 
perature, but saturates to a value easy to compute: at 
large T, (ie p goes to zero, the potential tends to a flat 
infinite well of size d, and thus the ground state reduces 
to ip(x) = ^2/dsm (xn/d) from which one finds imme- 
diately, for small x c , the saturation probability 
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There is therefore a minimum distance between necking 
of nanowires or localized collapses of nanotubes, given by 



2 l c d 3 

7T 3 X^X ' 



(13) 



and attained at very large temperature. This is a general 
feature of localized failures in Q1D systems whose lateral 
stability is protected by an infinite potential barrier. This 
is clearly not the case for DNA, where instead denatu- 
ration, or complete separation of the two filaments, is 
known to occur, as we will see later. No doubt this sat- 
uration effect could disappear should the defects come 
close enough to interact attractively-an interaction we 
completely neglect. Our estimates below show that this 
is not the case for carbon nanotubes. 

T c can be found by solving graphically for the eigen- 
value equation of the potential in Fig. [2j If A <C d, states 
with energy less than /3e p disappear at 



T c = A^/SA^/tt. 



(14) 



Numerical solution of the eigenvalue equation (Fig. [2]) 
shows that until about 0.8 T c no defects are present, yet 
1 — p grows steeply around T c [where it is roughly ( 1 — 
p(T c ))~(l-p) max /3]. 

While the phenomenology of nanowires is too rich to 



be considered here, we can use (13 1 for reasonable esti- 
mates in the case of carbon nanotubes collapse. If the 
radius is below a certain critical value (~ 3 nm), yet 
not too little (~ 1 nm), the cylindrical configuration is 
globally stable [5l [3] , yet the nanotube can collapse in a 
finite region. We first check our approximations. From 



( 14 1 we have that if A is less then an Angstrom, and 
k ~ e b ~ 1CT 1 eV A _x , then we obtain T c - 1CT 2 eV, 
or of the order of 10 2 -10 3 K. At very low temperatures 
the eigenstates reduce to plane waves in an infinite well 
of width A, and their eigenvalues are proportional to T, 
which is small: nevertheless, with the numbers above one 
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finds ei — eo ^> 1 if T 1 K, and thus the ground state 
approximation works at any practical temperature. Tak- 
ing x c = 0.34 nm (the equilibrium distance in the van 
der Waals interaction between graphene sheets), d = 2.8 
nm < d c (d c is the critical value above which the glob- 
ally collapsed configuration is stable, between 3 and 3.5 
nm [3 El E]) A ~ 1, l c ~ 1 nm, one finds the high tem- 
perature average length of non collapsed portions to be 



100 nm, or less. A few considerations: as d comes 



close to d c , the metastability of V(x) around x c becomes 
significant and the potential in Fig. [2] does not prop- 
erly describe the situation anymore. Conversely, as d 
gets smaller, Z m ; n also decreases, as smaller amplitude is 
needed in thermal fluctuations to induce failure, yet the 
potential barrier ef, increases, and with it the tempera- 
ture T c . Finally a reversed case exist for large (d > d c ), 
fully collapse nanotubes, which, conversely, we expect to 
show local thermally inflated regions. From ( 13 1 one can 
see that in this case, as d is larger, so is Z m i„, of the order 
of many hundreds of nanometers. 

Many theoretical and computational studies have been 
devoted to the study of bubbles in DNA [31 |5J d EH [ffl 
[T2~| . because of their importance for genetic transcrip- 
tion. In an inverted approach, now it is the joined double 
strand that we considered a "defect", and we ask our- 
selves what is the probability of existence of a bubble, or 
filament separation, of a certain length. We can use the 
formalism above by considering a square potential |14j . 
of width x c = 0.5 A and depth e& = 0.33 eV, [TU], roughly 
representing the chemical bonding between complemen- 
tary bases, on a semi-infinite line. As for the clastic con- 
stant we use k = 3 10 -3 eV/A 2 [TU]: that is consistent [8 



with the disappearance of a bond state and thus denatu- 
ration at a certain critical temperature T c = x c \/8keb/ir. 
With the numbers above, we can see that our single eigen- 
value approximation is justified for any reasonable tem- 
perature; therefore we find that the probability of a bub- 
ble decays exponentially with its length as Pyi — P , & 
result confirmed by computationally intensive numerical 
studies . Also, since p increases when the localization 
of the ground state in the potential well decreases, the 
average length of a bubble will be larger in the region 
where the bond between complementary bases is weaker, 
as expected for instance in regions rich of softer AT bases, 
a phenomenon demonstrated by recent work 18J. Also 
I increases with temperature, to reach the value p = 1, 
thus I — oo at the critical temperature for denatura- 
tion. For temperatures very close to denaturation, the 
single eigenvalue approximation might break down, and 
the probability might not then behave exponentially. 

In conclusion we have calculated the statistics of 
thermally- induced localized structural failures, defects or 
instabilities in Q1D systems. For the collapse of carbon 
nanotubes, and necking of nanowires, we predict a crit- 
ical temperature below which the occurrence of failure 
becomes extremely rare. We also show that the behavior 



of the average distance between defects is monotonically 
decreasing in temperature and approaches a minimum 
avearage distance for very high temperatures. These re- 
sults may be verified by experimental measurements of 
transport at different temperatures. For DNA, we predict 
that the probability of bubble openings decreases expo- 
nentially with the length of the bubble. It increases with 
temperature and is likelier to occur where the strength of 
bonding of complementary bases is weaker. These results 
may be verified by simulations using the PBD model. 

Discussions on DNA modeling with Boian Alexandrov 
(Los Alamos National Laboratory) was as useful as pleas- 
ant. This work was carried out under the auspices of the 
National Nuclear Security Administration of the U.S. De- 
partment of Energy at Los Alamos National Laboratory 
under Contract No. DE-AC52-06NA25396. 
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